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Abstract 

The evolution of the low-energy electromagnetic dipole response with the neutron excess is in- 
vestigated along the Sn isotopic chain within an approach incorporating Hartree-Fock-Bogoljubov 
(HFB) and multi-phonon Quasiparticle-Phonon-Model (QPM) theory. General aspects of the re- 
lationship of nuclear skins and dipole sum rules are discussed. Neutron and proton transition 
densities serve to identify the Pygmy Dipole Resonance (PDR) as a generic mode of excitation. 
The PDR is distinct from the GDR by its own characteristic pattern given by a mixture of isoscalar 
and isovector components. Results for the isotopes and the several N=82 isotones are 

presented. In the heavy Sn-isotopes the PDR excitations are closely related to the thickness of the 
neutron skin. Approaching ^'^'^Sn a gradual change from a neutron to a proton skin is found and 
the character of the PDR is changed correspondingly. A delicate balance between Coulomb and 
strong interaction effects is found. The fragmentation of the PDR strength in ^^^Sn is investigated 
by multi-phonon calculations. Recent measurements of the dipole response in ^^'^'^^^Sn are well 
reproduced. 

PACS numbers: 21.60.-n, 21.60.Jz, 21.10.-k, 23.20.-g, 27.40.+Z, 27.60. +j, 27.50.+e, 24.30.Cz 



I. INTRODUCTION 



The recent developments in experimental facilities of fast radioactive beams allows to 
study exotic nuclei far from the valley of beta-stability towards the neutron and proton 
driplines. One of the most interesting results was the discovery of a new dipole mode at 
low excitation energy. Typically, one observes in nuclei with a neutron excess, N > Z, a. 
concentration of electric dipole states at or close to the particle emission threshold. Since this 
bunching of 1~ states resembles spectral structures, otherwise known to indicate resonance 
phenomena, these states have been named Pygmy Dipole Resonance (PDR). However, only 
a tiny fraction of the Thomas-Reiche-Kuhn energy weighted dipole sum rule strength is 
found in the PDR region. Hence, these states will not alter significantly the conclusions 
about the importance of the Giant Dipole Resonance (GDR) in photonuclear reactions, as 
known for a long time 

Empirically, the low-energy PDR component increases with the charge asymmetry of 
the nucleus. The experimental situation from high-precision photon scattering experiments 
performed in neutron-rich stable nuclei has been reviewed recently in (2] • Of special interest 
are as well the newly performed experiments with radioactive beams in unstable oxygen 
3,Q] and tin isotopes, where observation of pygmy dipole strength has been reported Q. 

In jsl, 7, Q] we have established the PDR as a mode directly related to the size of the 
neutron skin. By theoretical reasons the PDR mode should appear also in nuclei with a 
proton excess [q]. Hence, the PDR phenomenon is closely related to the presence of an 
excess of either kind of nucleons. 

However, considering the measured dipole response functions the identification of a PDR 
mode is by no means unambiguous. While in a nucleus like ^ospb most of the dipole transition 
strength is found in the rather compact GDR region, the picture becomes more complicated 
in neutron-rich nuclei between the major shell closures. The dipole strength shows a ten- 
dency to fragment into two or several groups once the region of stable closed shell nuclei 
is left although the GDR still exhausts most of the (almost) model-independent Thomas- 
Reiche-Kuhn sum rule. A distinction between different dipole modes simply by inspection 
of the nuclear dipole spectra is clearly insufficient. As we will discuss in later sections a 
closer analysis of our theoretical results reveals a persistence of the typical GDR pattern 
of an out-of-phase oscillation of protons and neutrons in most of the low-energy satellites 
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until a sudden change happens in those parts of the spectra in the vicinity of the particle 
emission threshold. As we have pointed out before [g], Q, Q these states are special in the 
sense, that they are dominated by excitation of nucleons of one kind with small admixtures 
of the other kind. In a nucleus with neutron excess these excitations involve mainly neutron 
particle- hole (ph) configurations^. 

A clarification about the nature of the low-energy dipole strength in exotic nuclei can only 
be expected by the help of theory. With this paper we intend to contribute to this interesting 
problem by a systematic study of the evolution of dipole modes in the Sn isotopes. As 
suitable quantities, we consider the non-diagonal elements of the one-body density matrix. 
These transition densities give a snapshot of the motion of protons and neutrons during 
the process of an excitation. In principle, they are observables, e.g. for selected cases 
as the transition form factors in inelastic electron scattering (llj. In practice, however, 
experimental difficulties are usually inhibiting such measurements. 

Considering neutron-rich nuclei the character of the mean-field changes with increasing 
neutron excess, because of the enhancement of the isovector interactions. This has important 
consequences for the binding mechanism. In a neutron-rich nucleus the neutron excess leads 
to a rather deep effective proton potential, but produces a very shallow neutron mean- 
field. This results in deeply bound proton orbits with separation energies of the order 
of 20 MeV or even more as seen in empirical mass tables [12j]. The increase in proton 
binding is accompanied by a decrease in neutron binding. The most extreme cases are 

n n 

the spectacular halo states in light nuclei, e.g. |13l . 1141 ]. However, even for less extreme 
conditions unusual nuclear shapes are expected. In neutron-rich medium- and heavy-mass 
nuclei neutron skins of a size exceeding the proton distribution by up to about 1 fm have 
been predicted. Measurements in the Na p^] and Sn 16, [l3] regions indeed confirm such 



conjectures, although the present data are not yet fully extending into the regions of extreme 
asymmetry. 

The manuscript is organized as follows. In Section [III relations between a neutron or 
proton skins and the total nuclear dipole yield are discussed. In section IIIII the theoretical 
methods applied in this paper are explained. The mean-field part is treated microscopically 



[1] A certain amount of proton excitations is required as a compensation in order to suppress the motion of 
the nuclear center-of-mass. 
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by using HFB theory, but for the QPM calculations we allow empirical adjustments, which 
we incorporate by an phenomenological Density Functional Theory (DFT) approach. The 
description of response functions and transition densities by QRPA theory is discussed in 
section IIVI As an interesting mass region we explore the unstable tin isotopes and present 
results on the dipole response in section IVl Section |Vl] is devoted to comparisons to results 
of other calculations and data. The paper closes with a summary and an outlook in section 

m 



II. NUCLEAR SKINS AND THE DIPOLE RESPONSE 
A. The Dipole Excitations in Exotic Nuclei 

In this section we spend a few lines on discussing some relations between a proton or 
neutron skin and Pigmy Dipole Modes in exotic nuclei. We begin with recalling the definition 
of the nuclear electric dipole operator which in units of the electric charge and in terms of 

n 

the intrinsic coordinates is given by [18] 

^ = ^E^'^(i-^3.) = -^$^e>3. . (2.1) 

i i 

However, in nuclear structure calculations like the present one the particle coordinates 
are used which are related to the intrinsic coordinates by 



with the center-of-mass coordinate 



Be means of eq. 12.21 we find 



ir = ri-R (2.2) 



R = \Y.'^ . (2.3) 



^ = -^E(^3.-5)r-; (2.4) 

i 

= qnNRn + qpZRp (2.5) 

^Rp - Rn) (2.6) 



A 

where T3 is the 3-component of the total nuclear isospin operator which is a conserved 
quantity and can be replaced by its eigenvalue N — Z. The neutron and proton recoil 



corrected effective charges are denoted by g„ = —ZjA and = N/A, respectively. The 
partial sums over proton and neutron coordinates, normalized to the respective particle 
numbers, are denoted by -Rp,„. They describe the position of the center-of-mass of the 
protons and neutrons, but neither of the two quantities are conserved separately. 

As is obvious from eq. 12.11 only the protons participate actively in the radiation process 
while the neutrons follow their motion such that the position of the center of mass remains 
unperturbed as is reflected by eq. 12. 6[ This condition implies that R = ^Rn + ^Rp is a 
stationary operator. In other words, in any nucleus, whether stable or short-lived exotic, the 
condition R = must be fulfllled, leading to i?„ = —j^Rp. Together with eq. I2.6l this relation 
reflects the well known property of the dipole giant resonance (GDR) of an oscillation of the 
proton fluid against the neutron fluid, also found in hydrodynamical models of the nuclear 



dipole response 



19|. 



However, already in [19| it was pointed out that in a nucleus with a neutron excess 
Ne = A — 2Z > additional modes of excitation are possible. Such a system will also develop 
modes in which the excess nucleons are oscillating against the bulk, consisting of an equal 
amount of protons and neutrons, Zf, ~ A*";,. In that case we write R = ^Re + ^-^^Rb + ^Rp 
and another mode conserving the total center of mass is 

where Rb denotes the position of the center of mass of the remaining Nb = N — bulk 
neutrons. Because of Nb = Zb = Z we flnd the relation 

Re = —j^ (j^b + Rp^ , (2.8) 

indicating the motion of the excess neutrons against the core with an equal amount of 
neutrons and protons. Correspondingly, in a proton-rich nucleus, Zf, = A — 2N > 0, the 
excess protons may oscillate against the charge-symmetric bulk. In either case, watching 
that type of motion from the laboratory frame the core neutrons and protons will be seen 
to move in phase among themselves, but oscillate against the excess component, although 
the electric dipole operator D will couple directly only to the protons. 

If we assume harmonic motion, Rb,p,e{t) = -Rlpe^"*"^*' eq J2.8l leads to a relation among 



[2] We neglect higher order effects from the non-vanishing electric form factor of the neutron. 
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the amplitudes 




showing that in a neutron-rich nucleus the skin components will oscillate with an amplitude 
reduced by the factor compared to the core. 



B. Skin Thickness and Dipole Response 

The skin thickness is defined usually by the difference of proton and neutron root-mean- 
square (rms) radii 

6r=y/< rl>-^<rl>, (2.10) 

where 

.2 ^_ 1 / j3_2. 



< >=^ J d'rr'pgir) (2.11) 

denotes the rms radius of the proton and {q = p) and neutron (g = n) ground state density 
distributions pq, respectively, normalized to the corresponding particle number Ag = N, Z. 
For the present purpose a better suited choice is to express the differences in rms-radii in 
terms of the intrinsic coordinates, weighted by the 3-component of the isospin operator with 
eigenvalues ±1 for neutrons and protons, respectively, 

A,e = {o\j2^i^M = N<en>-z<ep> , (2.12) 

i 

which contains the same type of information as eq. I2.10[ 



If we express in terms of the laboratory coordinates r^, eq. 12.21 the ground state 
expectation value of eq. 12.121 leads to a corresponding expression in terms of the laboratory 
coordinates {r,}, 

A3r' = 5^<0|r3ir2|0> (2.13) 

i 

A . N-Z 2 



A-2 V A 
We have neglected contributions related to two-body correlations. 

In terms of the laboratory coordinates the intrinsic nuclear dipole transition operator 

could be expessed as 

D = Y,^i (gp^(l - T-3.) + g„^(l + raO) . (2.14) 



Defining the isoscalar (T = 0) and isovector (T = 1) charges and space vectors, respectively, 

qT = ^{qn + hfqp) ; XT = ^n(r3if (2.15) 

i 

we obtain the isospin representation 

D = qoXo + qiXi . (2.16) 

The reduced isovector/isoscalar dipole transition moments and the dipole transition proba- 
bihties can be expressed as 
=<0||(r3)^r1M>; 

B,{El) = \qoMP+q,MP\' . (2.17) 

For the present purpose we are interested in the isoscalar-isovector interference term, in 
particular 

^J2^d^ ■ ^d^* = (2-18) 

d 

\ d d d / 

Introducing a single particle basis {({){f)i\ enables us to express the nuclear dipole eigenstates 
\d> m terms of particle-hole excitations \a >. 

In the concrete case considered here, we use a QRPA description in terms of two- 
quasiparticle excitations \a >= \{ij)JM > with the single quasiparticle states i and j, 
respectively, coupled to total angular momentum JM. Neglecting ground state correlations, 
which are of at least 2p2h character, the left hand side of eq. 12. 181 can expressed in terms of 
single particle matrix elements 

lj:^p-M!p*=j:^s'-^t-^ (2-19) 

d i,j 

where vf = (|a^aj|) and = 1 — vf are occupation numbers. The last two terms in 
the above equation are of 2-body character and appear because of the pairing ground state 
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correlations. For our purpose the first term on the right hand side is of special interest. It 
is seen to correspond to 

■ M^f*v^ = J2<Af^J>-< Jlnf^t > (2.20) 
= < i\r- Tsr\i > = < i\r:i,r'^\i > v^, 

i i 

where we have used the completeness of the single particle states. Hence, we have derived 
an important theoretical relation between the non-energy weighted dipole sum rule and the 
skin measure defined before 

^■^^' = ^^{^B,(El)-4^\Mr? (2.21) 

d / i,j 

In the first term the pure isovector and isoscalar dipole sum rule strengths are subtracted 
off the full dipole sum rule, thus leaving the interference term. The isovector sum rule will 
be dominated, if not exhausted, by the GDR. From the Thomas-Reiche-Kuhn sum rule, i.e. 
the corresponding energy weighted dipole sum rule, we know that the GDR strength varies 
little along an isotopic chain, namely as NZ/A. With N = + 6N,A = Ah + 6N, and 
Zb = Z we find ^ ~ ^^^(^ + 'ifh^'^' emphasize again that the isoscalar sum rule is a 
pure recoil effect, appearing only in the neutron-rich nuclei and expressing the compensating 
motion of the neutrons in the laboratory frame. 

These relations reveal the intimate connection between the neutron skin (which is a 
static property) and the dipole spectrum (which is a dynamical property): From the above 
equations we find that the skin thickness is directly related to the dipole response. By eq. 
12.211 another aspect is emphasized, namely the fact that, seen from the laboratory, apparent 
isoscalar and isovector moments seem to contribute to the excitation of dipole states. As 
discussed above, modes involving an isoscalar component are allowed provided that the 
position R of the total nuclear center of mass is left untouched. 
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C. Dipole Excitations and Spurious States 



The problem posec 



well known, e.g. 



20 



_by broken symmetries in effective nuclear Hamiltonians like ours is 



2l| . For dipole excitations the most relevant effects are due to the 



violation of translational and Galilean symmetry. A known property of RPA is to restore 
the broken translational symmetry by generating a states at zero excitation energy, corre- 
sponding to a symmetry-restoring Goldstone-mode, provided that a complete configuration 
space was used. The transition strength scales with the total particle number A = N + Z 
and exhaust to a large extent the isoscalar sum rule. In practice, that is hardly achieved. 
But numerically we can enforce the restoration by a proper choice of the residual interac- 
tion in the isoscalar dipole channel. An alternative are projection techniques which, to our 
knowledge have never been applied to a realistic multi-configuration QRPA calculation. 

In this paper, we are considering the following situation. In a N ^ Z nucleus the condi- 
tions change insofar as new intrinsic excitations (see eq. fl2.8p ) will appear, not encountered 
in stable N c:^ Z nuclei. As already pointed out some time ago by Mohan et al. 191] in a 
neutron-rich nucleus the excess neutrons may be excited into oscillations against the core, 
either in phase or out of phase with the core protons. Especially the latter mode is the 
one from which we can expect a sizable content of isoscalar strength. That mode, however, 
will never appear as a pure isoscalar mode because of a compensating motion of the core 
neutrons required in order to keep fixed the center-of-mass of the whole system. Obviously, 
this is an intrinsic mode which will be strongly suppressed when approaching the N = Z 
line. In fact, the isoscalar content of the PDR states is impressively confirmed by a recent 
experiment in ^^°Ce {22], comparing spectra from inelastic scattering of a particles to (7,7') 
spectra. Since the a particle is a pure isospin T = probe it acts as an isospin filter and the 
spectra in 2^ show clearly the content of isoscalar transition strength in the PDR region. 



The special character of these transitions becomes clear from the shapes of the transition 
densities shown later which obviously do not resemble any expectations from classical or 
semi-classical models. 

The GDR mode as one of the most collective excitations in nuclei is well understood, both 
quantum mechanically and in semi-classical hydrodynamical approaches, while the nature 
of the PDR is still waiting for full clarification. The afore mentioned early attempts to 
incorporate the low-energy dipole modes into the hydrodynamical scenario [19] by a three- 
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fluid ansate seemed to work reasonably well in ^°^Pb, but when applied to the Ca-isotopes 23 1 
the model failed as pointed out by Chambers et al. [24| . Our more detailed microscopic QPM 
studies of the PDR strength in ^osp]^ including transition densities and currents, gave 
strong indications, that the PDR modes are of generic character, clearly distinguishable from 
the established interpretation of the GDR by strong vorticity components. The differences 
are also visible in the transition densities, where they are showing up in terms of a nodal 
structure, unknown from GDR excitations. Hence, we have the surprising situation that a 
mode with a more complex spatial pattern is seen at energies below the most collective state. 
This (theoretical) observation indicates that PDR and GDR states are indeed belonging to 
distinct parts of the nuclear spectrum. The situation is less confusing if we take the view 
that the PDR is related to a more complex excitation scenario as indicated by the transition 
densities and velocity fields discussed in [l^. The characteristic features of PDR transition 
densities will be investigated in the following for the whole chain of known Sn isotope, from 
io°Sn to i32Sn. 



III. PHENOMENOLOGICAL DENSITY FUNCTIONAL APPROACH FOR NU- 
CLEAR GROUND STATES 

A. The Density Functional 

Our method is based on a fully microscopic HFB description of the nuclear ground states 
and the quasiparticle spectra as the appropriate starting point for a single-phonon QRPA or 
multi-phonon QPM calculation of nuclear spectra. However, being aware of the deficiencies 
of existing density functionals, when leaving the region of stable nuclei we accept slight 
adjustments and phenomenologically motivated choices of parameters. We assure a good 
description of nuclear ground state properties by enforcing that measured separation energies 
and nuclear radii are reproduced as close as possible. 

We start by considering the ground state of an even-even nucleus in an independent 
quasiparticle model for which we use the microscopic HFB approach. The nucleons move in a 
static mean-field, which is generated self-consistently by their mutual interactions including 
a monopole pairing interaction in the particle-particle (pp) channel. Following the DME 



approach of 



261], the interactions are taken from a G-Matrix, but renormalized such 
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that nuclear matter properties are reproduced, thus accounting effectively for correlations 
missed by a static two-body interaction. In local density approximation the problem is 
then reduced essentially to the level of a Skyrme-HFB calculation as discussed in 25||. An 
important difference, however, is the use of a microscopically obtained density dependent 
pairing interaction. The HFB and BCS equations are solved self-consistently with state 
dependent gaps by iteration, until convergence of the mean-field and single-particle energies, 
gaps and densities is achieved. 

The single-particle energies and ground state properties in general are critical quantities 
for extrapolations of QRPA and QPM calculations into unknown mass regions. Here, we 
put special emphasis on a reliable description of the mean-field part, reproducing as close as 
possible the g.s. properties of nuclei along an isotopic chain. This is achieved by solving the 
ground state problem in a semi-microscopic approach. Following the arguments given in 6| 
we take advantage of the Hohenberg-Kohn 27|] and Kohn-Sham 28|] theorems, respectively, 
of density functional theory, which state, that the total binding energy B{A) can always be 
expressed as an integral over an energy density functional with a (quantal) kinetic energy 
density r and density dependent self-energy parts U{p), respectively, 

B{A) =Y. j U{p) + \p,U,{p)\ + , (3.1) 

q=p,n ^ ' 

where we have chosen a representation in terms of proton (g = p) and neutron (g = n) 
densities pq = Pq{r), respectively, as appropriate for nuclei far from the stability region with 
exotic charge-to-mass ratios. The total isoscalar (T = 0) and isovector (T = 1) densities are 
defined by pr = Pn + {—)^Pp and po is normalized to the total particle number A. 

In addition to the kinetic and potential energy terms eq. ( 13. ip includes pairing contribu- 
tions, which are indicated separately by E^""^^ . In fact, this means, that we use an extended 
version of the Kohn-Sham theorem including the proton and neutron pairing densities Kq 
as well, as dictated by HFB theory. Hence, the density functional underlying in eq. (13. ip 
is of the form £{t, p, k), where each of the arguments are understood to include proton and 
neutron parts, respectively. 

In terms of the single-particle wave functions ^Pjq{r) and the occupancies Vj^ the kinetic 
energy density is given by 
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and the number and pairing densities are 

p^(f) = 5^t;|>,,(f)p (3.3) 



l^VjgUjg\ifijg{r)\'^, (3.4) 



2 
J 

where vjg, ujg denote BCS amphtudes with u^^ = 1 — fj^. The summations over j includes 
the full set of quantum numbers specifying the single-particle states ipjgif). 

Rather than using a conventional density functional like the Skyrme functional we choose 
to express the interaction part in terms of a superposition of central and spin-orbit potentials 
of Wood-Saxon shape. This ansatz gives us the full flexibility to describe nuclear ground 
state properties like binding energies, root mean square radii, and separation energies to the 
required accuracy. The price to be paid is a lack of contact to a fully microscopic picture like 
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251]. However, for the present purpose and in view of the persisting uncertainties on the 
dynamics in strongly asymmetric nuclear matter, we are convinced, that a phenomenological 
approach allowing a self-consistent description of nuclear ground states is an eligible method. 

Hence, in order to describe the bulk properties of the nuclear ground states in the best 
possible manner, we decide to be satisfied by using functionals optimized to a given mass 
region, in this case the Sn isotopes. The parameters of the model are the strengths, the radii 
and the diffuseness parameters of the corresponding parameters. A posteriori the collected 
information will allow us to derive eventually a nuclear energy density functional of general 
applicability. In other words, we try to avoid a biased choice of operators by assuming a 
certain operator structure at the level of a two-body interaction. 

B. The Single-Particle States 

From eq. ( 13. ID we derive by variation a Schroedinger equation 

(-^V^ + S,(r-') - r^.g^ y,(f) = (3.5) 

for the single-particle wave functions (pjg and eigen energies rjjg. The self-energy appearing 
in eq. (13.51) is obtained variationally from the interaction energy density 

Eint = ^^PqUg{p), (3.6) 
Q 
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where we have defined the single-particle occupation probabilities p^, which in BCS approx- 
imation are given by f^. By variation with respect to pq we obtain 

S.(P) = ^^ Y.P'i'U'^'^p) ■ (3-7) 

Because of the intrinsic density dependence of Uq{p) we find, that differs from the proper 
interaction energy by a rearrangement potential 

^,{P) = UM + U^;\P) (3.8) 

given by 

which is discussed in more detail in App. |X1 In nuclei with non-vanishing pairing addi- 
tional contributions from the density gradients of E^^"'^'^^ = E^°'^^\n^ p) will also contribute. 



C. Pairing and Quasiparticle States 

From the density functional, eq. (13.11) . we obtain the proton and neutron pairing fields 
Aq(p, k) by variation with respect to the pairing (or anomalous) densities Kg, eq. (13. 4p 

^{p,K) = ^-^ = KqV^^^'^\p), (3.10) 
dKq 

which we decide to factorize into the anomalous density and a local density dependent 
pairing strength V^*^^"*'')(p), depending on the local bulk density p = p{r). V^^^-^^^p) is 
discussed below. 

With the usual Bogolubov transformations we obtain the quasiparticle states = 



29|]. Together with the Schroedinger equation (13. 5p . we solve self-consistently 
the BCS gap equation for the state dependent pairing gaps Ajq = {jq\Aq\jq) for protons 
and neutrons, respectively, 

A,, = J2'^kqU,'^(f^\V^'^'^^\k~k) = . (3.11) 

k ^ k ^^'i 

Time-reversed states are denoted by a tilde. In a spherical symmetric nucleus the BCS state 

amplitudes are 



2 V E,q 
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with the quasiparticle energy 



En=\Km,-\Y + ^l (3.13) 



The pure mean-field single-particle energies rjjq are obtained from eq. (13.51) . The proton and 
neutron chemical potentials are denoted by Xp ^^ respectively. 

In the practical HFB calculation we use a pairing strength of a simple form 

V(p-^)(p) = (v,^,[l - (^)/5) + VU-f^ , (3.14) 

simulating the on-shell singlet-even NN interaction amplitude and depending on the local 
nuclear density p = p{r). The interaction strength Vext = —9280MeV fm^ is determined 
such, that asymptotically for p — > the nucleon-nucleon scattering length app ~ ann ~ (Ise = 
— 17.3/m in the singlet even channel for like particles is reproduced. Vint = —0.721MeVfm^ 
is fixed by requiring, that the pairing gap has a maximum of A(pc) = 2 MeV at = l/3po of 
the equilibrium density of infinite nuclear matter, po = 0.16/m^'^. The best results are ob- 
tained for a small value of the density exponent = ~ 0.012. In the pairing calculations 
we include proton and neutron single particle states up to the respective continuum thresh- 
olds. In this way, we avoid instabilities in the BCS equations and the calculations of number 
and pairing densities due to the possible admixture of unbound quasiparticle orbitals into 
the bound state region. Such an approach is permissible because in all the consider nuclei 
the driplines are not reached. Hence, the more involved treatment by explicitly solving the 



coupled Gorkov-equations, discussed e.g. in [30|, can be avoided without a significant loss 
of accuracy. 



D. HFB Results for Sn Isotopes 

As discussed above, we decide to express the full proton and neutron self-energies Sp „ in 
terms of (a superposition of) Wood-Saxon potentials Spl^"^-* by a least-square fit of the depth, 
radius and diffusivity parameters to separation energies and charge radii, taken either - if 
available - from empirical mass compilations [12!], or from our HFB calculations. Different 
to the usual HFB approach the s.p. wave equations are solved with effective mass m*=m 
thus removing the known problem of unrealistically large HFB level spacings at the Fermi 
surface. 
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FIG. 1: (Color online) Ground state properties of the Sn isotopes. The nuclear binding energies 
per particle calculated with the DFT approach, discussed in the text, are compared to data from 
the Audi-Wapstra compilation [l2]. 

The reproduction of the total binding energy B{N,Z), calculated as indicated above, 
of the charge radius and the (relative) differences of proton anc 



neutron root-mean-square 



(RMS) radii 5r, taken from our previous HFB calculations |25U26l|. are imposed as additional 
constraints. The results of the represent approach are displayed and compared to measured 
values in Fig. [H The ground state neutron and proton densities are displyed in Figl2] for 
several tin isotopes. The comparison between the neutron and proton densities, obtained by 



HF calculations with the D3Y G-matrix interaction (see Fig. 7 from ref. 25|]) and the present 
ones is very reasonable. 

Of special importance for our investigation are the surface regions, where the formation of 
a skin takes place as is visible in FigI21 For A>106 the neutron distributions begin to extend 
beyond the proton density and the effect continues to increase with the neutron excess, up 
to ^^^Sn. Thus, these nuclei have a neutron skin. The situation reverses in ^^^^^^"^Sn, where 
a tiny proton skin appears. 
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Radial distance [fm] 



FIG. 2: (Color online) BCS ground state densities of Sn isotopes obtained by the phenomenological 
DFT approach and used in the QPM calculations. 

IV. QPM DESCRIPTION OF NUCLEAR EXCITATIONS 
A. The QPM Hamiltonian 

The excitations are calculated in the framework of the Quasiparticle-Phonon Model 



(QPM) with the model Hamiltonian 291]: 



H = HMF + Ht; + Hf^ + HZ . (4.1) 

Here, Hmf = Hgp + Hpair is the mean-field part, discussed in the previous section. Hence, 
different from the standard QPM scheme we use single-particle energies and wave functions, 
obtained self-consistently, according to the procedure described above. For the QPM cal- 

16 



culations the pairing part is simplified by using a constant matrix element. The method 
we have applied for the determination of the ground state properties has been successfully 
applied fo. the investigation ot low-lyn.g d.pole ntodes in the tn. .otopic cha.n before fiQ 
and more recently also in the = 82 isotones |8| . 

H^, H^sl/i ^'^'^ -^M residual interactions, taken as a sum of isoscalar and isovector 
separable multipole and spin-multipole interactions in the particle-hole and multipole pairing 
interaction in the particle-particle channels. The latter is included only for the quadrupole 
and octupole excitations. 

Building blocks of the model space are the Quasiparticle-Random-Phase-Approximation 
(QRPA) phonons: 



hi2 

defined as a linear combination of two-quasiparticle creation and annihilation operators 
A\fj_, respectively. The latter is the time reversed operator = {—)^^^A\^fj_. 

Here j = [nljmT) is a single-particle proton or neutron state. 

The (bare) two-quasiparticle operators 



are defined by coupling the one-quasiparticle operators to total angular momentum A with 
projection fj, 




(4.2) 




(4.3) 




(4.4) 




{jimij2m2\Xfi). 



mxMJO) = 1, 



(4.5) 



which can be rewritten in terms of two-quasiparticle weight factors 




(4.6) 



jl>j2 



iij2 1 



iii2 1 
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The weight factors Wj^j^iX^i) are given for some states in Table I and Table II, respec- 
tively. 

The QRPA operators obey the equation of motion 

[H,Qt]=E^Qt , (4.7) 

which solves the eigenvalue problem, giving the excitation energies Ea and the time-forward 
and time-backward amplitudes [29i] i^'jlj^ and '^^Ij^-, respectively. 

The QPM Hamiltonian (14. ip is rewritten in terms of phonons {2^: 

H = Hph + Hgph = ^ ExiQX^iQxi^i (4.8) 



+ 9 XI C'AimA2M2^Aiir(^3^3) 



2 

AiA2A3Jil2i3/ll/i2M3 

The first term in the eg. (14.81) contains free phonon operators and refers to the harmonic 
part of nuclear vibrations, while the second one is accounting for the interaction between 
quasiparticles and phonons. The latter reflect in anharmonic effects and fragmentation of 
the nuclear excitations. 

The Hamiltonian (14.81) is diagonalized assuming a spherical 0^ ground state which leads 
to an orthonormal set of wave functions with good total angular momentum JM. For even- 
even nuclei these wave functions are a mixture of one-, two- and three-phonon components 



3l| in the following way: 

+ E (^^) [Qt..n X ^A2.2 J JM + E 



(4.9) 



JM 



*0 



Alii x-^i 
A2i2 A 

2.AiiiA2i2//j A rrg+ .(g)Q+ .1 . 

A313 V-" J Ai/iili ^ A2/i2«2J ^ ^ A3/i3l3 

where R,P and T are unknown amplitudes, and v labels the number of the excited states. 

The nuclear response on an external electromagnetic field is described in terms of quasi- 
particles and phonons by a transition operator composed of two parts: 

M(E(M)A/x) = MP'^(E(M)A/i) + M''P'^(E(M)A/i) (4.10) 
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The first part is responsible for tlie transitions witli one-plionon exchange between the initial 
and final states. The second one contains structures [a^ ® including the interaction 

between quasiparticles and phonons. It is important for the description of the so-called boson 
forbidden transitions between nuclear states with the same number of phonons or differing 
by an even number of them. The corresponding equations of each of the terms could be 
found in [32|. 



B. Transition Densities 

In order to understand the character of a nuclear excitation it is useful to consider the 
spatial structure of the transition. This is accomplished by analyzing the one-body transi- 
tion densities 6p{r), which are the non-diagonal elements of the nuclear one-body density 
matrix. Physically, 5p{r) corresponds to the density fluctuations, induced by the action of 
an (external) one-body operator on the nucleus. Hence, the transition densities are directly 
related to the nuclear response functions and by analyzing their spatial pattern we obtain a 
very detailed picture of e.g. the radial distribution and localization of the excitation process. 



io|. 



The particular usefulness of such an analysis for PDR states was pointed out in 

Using the complete set of single-particle states ^j{r) from Hmf and a multipole expansion 
by means of the Wigner-Eckardt theorem, we find the isoscalar (T=0) and isovector (T=l) 
transition densities in second quantization: 

5p^{r)= J2 (^'>^A.(r))Vj;,(r) [a>,4^ . (4.11) 
iii2;AAi 

For the present purpose we consider non-spin flip transitions of isoscalar and isovector 
character. The radial parts are given by binomials of radial single-particle wave functions 
and reduced matrix elements 

P?3.W = ^;i(0%(0i(jill^'n||j2)(g|r3^|g) , (4.12) 

A 



with A = v2A+T. The isospin matrix element {q\Tr[\q) is unity for T=0. For an isovector 
transition we have {q\T^\(l) = ±1 for neutrons and protons, respectively. 

The transition densities are obtained by the matrix elements between the ground state 
|\l>i) = \JiMi) and the excited states = \JfMf), 

pU^= E (^'^A,(r))V^S.W(^.^^^/lK«..]Aj'^^^^) (4-13) 
iii2;AM 
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Here, we are interested only in the two-quasiparticle creation and annihilation parts which 
are given by 

r^^(jij2) = {u,,Vj, + v,,u,,) (a+^(jiJ2) + A^^ijm)) (4.14) 
Equation (I4.14p can be rewritten in terms of QRPA phonons defined by the relation (14. 2p : 

rtOU2) = Y.3nn {QU + ^A,.), (4.15) 

i 

where 

9nn = ^71^^^^^^'^ iVn + Vn) (4-16) 
accounts for the BCS and QRPA properties, respectively. Thus, we find 

pU^= E (^'^A,(r))V-3.(r)(J/M,|r+^(jij2)|J.M,) . (4.17) 
iii2;AM 

We identify iJjMj) = |0) with phonon vacuum and obtain the excited states by means of 
the QRPA state operator, eq. (14. 2p . \JfMf) = Q^^JO), which leads us to the commutator 
relation 

pU^= E (^'>^AM(r))Vj;,(r)(0|[g,,„r+^(jij2)]|0) . (4.18) 

Hence, in QRPA theory the one-phonon transition density is given by the coherent sum over 
two-quasiparticle transition densities entering in the structure of a phonon by the relation: 

The shape of the transition density defined by eq. (14.190 is rather strongly correlated with 
the collectivity of the phonon. For example, the transition densities of the non-collective, 
two-quasiparticle excitations typically have pronounced maxima inside the nucleus. Those 
corresponding to the collective transitions with a large number of coherently contributing 
two-quasiparticle transitions have a maximum at the nuclear surface. 

The reduced transition probability B{EX) for the excitation of a state Jf from the ground 
state Ji is connected with the transition density with the relation: 

2 



BiEX) 



LT=0 



1 /"OO 



(4.20) 



2J, + 1 

where denotes the effective isoscalar and isovector charges, respectively, introduced before. 



20 



C. The QPM Model Parameters 



Following ref. [33|, the ratio k[^^ /kq^^ of the isovector and isoscalar multipole strength 
parameters, respectively, is assumed to be a constant, independent of the multipolarity A. 
We can find this ratio from the dipole coupling constants by projecting the spurious 1~ 
state to zero excitation energy and fitting the experimental energy of the Giant Dipole 
Resonance (GDR)^, [s^. For those nuclei, where GDR data are not available the empirical 
^max^ = law is used. The El transition matrix elements are calculated with recoil- 

corrected effective charges g„ = —Z/A for neutrons and qp = N/A for protons, respectively, 
as discussed in sect. HTl 



V. RESULTS FOR THE DIPOLE RESPONSE 
A. General Features of the Dipole Response 

The dependence of the calculated total PDR strength B{E1) f) on the mass number 
for the whole chain of isotopes ^°°~^^^Sn is shown in Fig 13] and compared to the skin thickness 
6r, eq. (12.101) of these nuclides. Here the sum is taken over QRPA dipole states presented 
in Table 1 and Table 2, respectively. According to state vectors structure they have been 
associated with PDR. These results illustrate and confirm the conclusion drawn in sect. 
ITT] and as already stated in previous work Q establish the close relationship of the PDR 
strength and the skin thickness. 

In the region between ^^°~^^^Sn the total PDR strength increases smoothly with the 
neutron number. This establishes a clear correlation of the total PDR strength and the 
thickness of the neutron skin in these nuclei, thus confirming our previous results 
and the more recent investigations for several = 82 isotones {s]. The close relationship 
between 6r and the PDR modes is underlined by the result, that the PDR becomes negligibly 
small in the region ^°^^^°^Sn, where 6r changes sign (see FigjS]). In these isotopes the lowest- 
lying states carry the characteristics of the low-energy branch of the GDR as indicated by 
the structure and shape of the transition densities. 



[3] Note, that in our previous work 6] the dipole response in ^^'^ ^"^^Sn was calculated with the bare charges, 
leading to systematically smaller values of the total transition strength. 
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FIG. 3: (Color online) QPM results for the total PDR strengths in the ^^^Sn isotopes (upper 
panel) are displayed for comparison together with the nuclear skin thickness 5r, eg. (j2.10p (lower 



36( and ^^^s 



panel). Experimental data on the total PDR strengths in ^^^Sn and ^-^^Sn ref. 
are also shown. In the lower panel, the skin thickness derived from charge exchange reactions by 
Krasznahorkay et al. Q, are indicated. 



From the QRPA calculations in ^^°^^^^Sn a sequence of low-lying one-phonon 1~ states 
at excitation energies E* = Q — 7.5 MeV of almost pure neutron structure is obtained with 
a minor fraction of protons less than 1%. The structure of the state vectors is indicated 
in Table I. The most important part of the total PDR strength comes from the excitations 
of the least bound neutrons from the 3s- and 2p- and 2d-subshells. Some other neutron 
orbitals of significance for the size of the neutron skin are l/iii/2 and lii3/2, which have an 
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important contribution to the PDR transition matrix elements in ^^^Sn nuclei, for example. 

The dominant neutron structure and remarkable stability of the wave functions of the 
low-lying one-phonon 1~ states in these nuclei is in agreement with our previous findings on 
the PDR mode in ^^'^"^'^^Sn isotopes 6|] and N=82 [8] isotones. 

Towards the lighter Sn isotopes the average energy of the excited dipole states increases, 
while their total number decreases (see Table I). The dependence of the PDR energy on the 
neutron excess is connected with the one-neutron separation energy, decreasing gradually 
towards the heavier tin isotopes. Such a tendency has been observed also experimentally in 
N=82 isotones (see ref. js^). 

An exception is the double magic ^^^Sn nucleus, in which the PDR centroid energy is 
E* = 7.1MeV and is still below the neutron emission threshold. The present result in 

Sn is slightly different from our previous one due to minor readjustments in the 

single-particle spectrum. 

In ^''^"^''^Sn the lowest dipole excitations, E*= 8.1-8.3 MeV, are dominated by proton 
excitations. The structure of the QRPA state vectors and B(E1) transition probability 
are given in Table II. There it is seen, that configurations involving quasibound 2j93/2 and 
1(79/2 proton states confined by the Coulomb barrier are the major components. This is 
a remarkable Coulomb effect enlightening the delicate balance among various effects as a 
prerequisite for a PDR and, by comparison to Fig. [2], the existence of a nuclear skin. From 
Fig. [2] it is seen, that this is the mass region, where the neutron skin turns into a proton 
skin. In agreement with the considerations in sect. [Ill the vanishing skin is accompanied by 
a strong suppression of the dipole strength. The smallest strength is found at A = 110 — 112, 
which is slightly above the turnover point of 6r at A = 106. This delay is caused by Coulomb 
effects, which enhance the dipole response from weakly bound proton orbitals in that mass 
region over the values to be expected for full isospin symmetry. 

Electromagnetic breaking of isospin symmetry is also the main reason for the persisting 
of low-energy dipole strength close to ^^''Sn. Already the quite different behavior is an indi- 
cation for another mechanism underlying these excitations. There, at N = Z the isoscalar 
dipole charge vanishes, hence the electromagnetic operator by itself does no longer support 
isoscalar transition. However, Coulomb effects in the single-particle wave functions translate 
into an intrinsic isospin symmetry breaking on the level of matrix elements. The mechanism 
behind a neutron skins in the heavy Sn isotopes is a strong interaction effect, namely the 
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repulsive action of the isovector self-energy to the neutron mean-field. In neutron-rich nuclei 
the isovector self-energy adds attractively to the proton potential, which partially compen- 
sates the Coulomb repulsion. Because for N— the isovector self-energy becomes negligible 
the proton skins seen in the light Sn isotopes must be of a different origin. In fact, they 
are due to the Coulomb potential. Towards N = Z the Coulomb interaction can act in full 
strength on the protons, pushing them apart and leading to a rearrangement of a certain 
fraction of the nuclear charge in the surface region. 

The average energies in Table I and Table II have been obtained by the relation 
{E) = JZi^i^i/ J^i^ij where Ei and Bi are the QRPA energies and reduced transition 
probabilities, respectively. 

B. PDR and GDR Transition Densities 

For a more detailed insight into the characteristic features of the dipole excitations we 
consider the evolution of the proton and neutron transition densities for in the various energy 
regions. In FigH] and Figl5] we display the QRPA transition densities for several = 82 
isotones in and the Z = 50 isotopes 112,122, i32gj^ three different regions of excitation 
energies: the low-energy PDR region below the neutron emission threshold, the transitional 
region up to the GDR and in the GDR region and beyond^. The transition densities 
displayed in Fig HI Fig. [5] and Fig. [6] were obtained by summing over the transition densities 
of the individual one-phonon states located in the energy intervals denoted at the top of 
each column of the figures, i.e. 



Pxi{r) is determined by equation fl4.19p . where the module and the phases are unambiguously 
determined by our microscopic approach. The neutron and proton transition densities are 
then obtained by taking half the difference and the sum of the isoscalar and isovector pieces, 
respectively. 

A common features of the all cases presented in Fig.4 is that up to E* = 8.1 MeV the 
protons and neutrons oscillate in phase in the nuclear interior, while at the surface only 
neutron transitions contribute. The same behavior of the neutron and proton transition 

[4] For a detailed discussion of the dipole response of the N=82 isotones we refer to ref. 




(5.1) 
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E*= 0-8.1 MeV E*= 8.1-8.6 MeV E*> 8.1 MeV 
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FIG. 4: (Color online) QRPA results for the one-phonon dipole transition densities in N=82 nuclei. 

densities is observed below 8 MeV for tin isotopes (see Fig.5). This pattern is generic to 
the lowest dipole states making it meaningful to distinguish these excitations from the well 
known GDR states. Hence, we are allowed to identify the PDR states with a new mode of 
nuclear excitation, not seen in stable N ^ Z nuclei. 

In the energy region E* = 8.1 — 8.6 MeV for N=82 isotones and 8-8.5 MeV for Z=50 re- 
spectively, the transition densities suddenly change. Rather abruptly, protons and neutrons 
start to oscillate out of phase over the whole nuclear volume as known from the GDR. Thus, 
we are encountering the low-energy part of the GDR, although the strengths of these two 
different type of excitations, the PDR and the low-energy GDR tail are quite comparable. 
Also energetically they are located very close to each other. This makes the task to dis- 
tinguish the two modes rather demanding. Theoretically, we can always use the transition 
densities for a detailed analysis and a precise identification of the mode although a corre- 
sponding experimental measurement will not be feasible for the foreseeable future. Finally, 
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FIG. 5: (Color online) QRPA results for the one-phonon dipole transition densities in Ii2,i22,i32gj^ 
nuclei. 




FIG. 6: (Color online) QRPA results for the one-phonon dipole transition densities in ^'^"Sn. 
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we show in FigH] and Figl5] the one-phonon QRPA proton and neutron transition densities 
for the states at E* > 8.1 MeV and E* > 8 MeV, respectively. A pronounced isovector 
oscillation of protons against neutrons, peculiar to the GDR, is observed. The latter, as a 
very collective mode, has a strength one order of magnitude larger than the PDR. 

We have pointed out the special character of the low-energy dipole excitations in the 
ioo-io4gj-^ isotopes. This is reflected also in the transition densities. In ^°°Sn, FiglHl proton 
oscillations prevail below E* < 8A MeV. In the nuclear interior isoscalar, or mixed symmetry 
vibrations of protons and neutrons are found, while at the surface only protons contribute. 
Hence, this mode could be related to a proton skin excitations. In the energy region E* = 
8.4 — 9.2 MeV oscillations of weakly bound neutrons from the surface region take place. The 
behaviour of the transition densities and the structure of the 1~ states at these energies is 
similar to the neutron PDR mode identified in the more neutron rich tin isotopes (see Figl5] 
as well). At energies E* > 9.2 MeV the low-energy tail of the GDR is encountered. The last 
plot in FiglH] displays the neutron and proton transition densities summed over the GDR 
region, E* = 9.2 - 20 MeV. 

The QRPA calculations on the neutron and proton PDR and the GDR strength distri- 
butions at excitation energies up to 20 MeV in several thin isotopes in the mass region 
^°''Sn-^^^^Sn are presented in Fig. ([7]). An interesting feature we have observed is that be- 
tween the N=50 and N=82 closed shells with the increase of the neutron number from ^°°Sn 
toward ^^^Sn the PDR strength is shifted to lower excitation energies relatively to the GDR 
mode, which is almost unchanged. This can be explained with a strong correlation between 
the PDR excitation energy and the energy of the neutron threshold, which also decreases in 
the same direction. 



C. Multi-phonon effects in the low-energy dipole spectra of Sn. 

In the multi-phonon QPM calculations the structure of the excited states is described 
by wave functions as defined in eq. l4.9[ We now investigate multi-phonon effects using a 
model space with up to three-phonon components, built from a basis of QRPA states with 
J'^ = 1^, 2"*", 3~, 4+, 5^, 6+, 7^, 8+. Since the one-phonon configurations up to E*=20 MeV 
are considered the core polarization contributions to the transitions of the low-lying 1~ states 
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TABLE I: Energy, B(E1) values and wave functions of the first QRPA 1~ states in the "^^'^""^^^Sn 
isotopes. Only the dominant neutron and proton components are given. Neutron and proton 
configurations are denoted by the indices v and tt, respectively. 



Nucleus 


State 


Energy 


Structure 


B(El)t 


<E> 




K 


[MeV] 




pfm^l 

L J 


[MeV] 

L J 


iiOSn 


ir 

i 


7.834 


99.9%[l07/22/7/2]z/ 


0.001 


7.8 


ii^Sn 


ir 

i 


7.509 


99.8%[l07/22/7/2lz/ 


0.001 








7.906 


99.0% [3si/23p3/2l 


0.144 


7.9 


kJll 




7 '^99 


Q9 8%f1 n., ;o9 f-, /0I7/ 


001 








7.665 


99.2% [3si/23p3/2l V 


0.159 


7.7 






8.021 


99.9% [2d3/23!?3/2l 


0.005 




ii^Sn 


ir 


6.974 


99.7%[207/23f7/2lz^ 


0.001 






1 - 

^2 


7 188 


yy./o [•J''i/2'J/'3/2j z^ 


1Q9 


7 9 








+ 0.1% [laortl/illrtl TT 








1^ 

o 


7.391 


99.9% [2(i3/23pi/2l 


0.009 




iiSSn 


ir 
i 


6.904 


99.6%[l07/22f7/2lz/ 


0.001 






1 ~ 

^2 


7 054 


Q8 1% r^Si /o^no/ol 7y 
»o.-L/o [»>''l/2*'/'3/2j 




7 1 








+ 0.1% ri(7q/2l/ll1/2l TT 








1^ 

o 


7.098 


99.1% [2(i3/23p3/2l 


0.02 




i20Sn 


ir 

i 


6.795 


99.6%[2o?3/23p3/2li/ 


0.009 








6.870 


95.2% [197/22/7/2] 


0.014 


6.9 




I3 




6.910 


94% [3Si/23p3/2] 


0.238 










+ 0.1% [159/21/^11/2] ^ 






i22Sn 




6.469 


99.8%[2d3/23p3/2]z/ 


0.014 








6.710 


95.3% [3si/23p3/2] u 


0.245 


6.7 








+ 0.1% [I59/2I/III/2] TT 










6.754 


95.8% [157/22/7/2] 


0.009 










+ 0.1% [139/21/^11/2] TT 
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Table I continued 








ir 


6.359 


99.8%[2d3/23p3/2]j^ 


0.017 








6.702 


94.8% [3si/23p3/2] V 


0.284 


6.68 








+ 0.1% [159/21^11/2] TT 








I3 


6.749 


95.6% [1^7/22/7/2] 1^ 


0.009 










+ 0.1% [lff9/2l/ill/2] ^ 






Sn 




6.180 


99.7%[2(i3/23p3/2]i^ 


0.019 






^2 


6.621 


51.4% [3si/23p3/2] ly 


0.163 


6.6 








+ 48.3% [I57/22/7/2] 








la- 


6.642 


51.5% [1^7/22/7/2] 1^ 
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+ 47.% [3si/23p3/2] 1^ 
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ir 


5.611 


99.7%[2d3/23p3/2]i/ 
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6.201 


97.8% [3si/23p3/2] 1^ 


0.306 


6.2 
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6.352 


99.1% [I57/22/7/2] 


0.001 




Sn 


ir 
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^ 
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98.1% [3si/23p3/2] ly 
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6.114 


99.4% [I57/22/7/2] 
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5.754 
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88.6% [24/22/7/2] ly 
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+ 0.2% [159/21/111/2] TT 







are taken into account explicitly. Hence, we do not need to introduce dynamical effective 
charges. In the excitation energy interval up to E*=9 MeV we use a total of about 250 
multi-phonon configurations. 
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TABLE II: The same as Table I for i°°^i°''Sn isotopes. 



Nil 1 r*! pn t; 


State 


i^jiiei 


S^f Vl 1 1 1 TP 




< F, '> 




TTT 


[MeV] 


w ■ % 


[e^fm^] 


[MeV] 


loogn 




8.032 


99.5%[l/5/22(i5/2]l/ 


0.001 








8.292 


82.1% [2p3/22(i5/2] vr 


0.028 


8.29 






8.174 


82.6%[2p3/2 2d5/2]7r 


0.031 








8.256 


80.8%[2p3/22d5/2]7r 


0.016 





The results for the dipole response below the neutron threshold in ^^''Sn are presented 
in FiglHl By comparing FiglH^ and FigJH^ it is seen that the pure two-quasiparticle QRPA 
strengths is strongly fragmented once the coupling to multi-phonon configurations is allowed. 
As found previously the lowest-lying 1~ state without a QRPA counterpart is predominantly 
given by a two-phonon quadrupole-octupole excitation [2,]. The [2^'" (g) 3|f] configuration 
accounts for 85% of the QPM wave function. The 1^ state is located at E* = 3.50 MeV, 
carrying a reduced transition probability B{El;g.s —>■ li) = 6.06 10"'^ e^fm^. The values 
are in a good agreement with the experiment {E* = 3.49 MeV and B{E1; g.s — ^ 1^^) = 6.08 



3.0,0. 



10^^ fm? [39] and previous QPM calculations refs. 

Here, our attention is especially focused on the 1~ states above the two-phonon dipole 
state and below the neutron threshold. From the analysis of the QRPA calculations discussed 
above the 1~ states presented in FigJH^ are PDR modes. Their fragmentation over the 
multi-phonon 1~ excited states are shown in FiglHjD. From that plot we can determine the 
energy region, where the PDR is located. In the particular case of ^^^Sn in the interval 
E* < 7.5 MeV is exhausted about 80% of the total one-phonon PDR strength. 

Comparing the fragmentation pattern of the theoretical low-energy dipole strength in 



^^^Sn to recent measurements 
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391], displayed in FiglHli, we find that the three-phonon 



QPM results are still not fully accounting for the observed distribution, although we have 
increased the number of^he multi-phonon configurations twice in comparison with our pre- 
vious calculations from [6]. However, the calculated total QPM dipole strength in the PDR 
energy range E* = 5.7 — 7.2 MeV is B{E1)qpm = 0.324 e^fm? which is almost identical 
to the experimentally deduced strength, Yl, Bexp{,El) = 0.345(43)e^/m^. 



30 



10 


1 


1 


1 


10"' 


1 


1 U 




1 u 




10 


-j 


1 


-! 


10-' 


i 


10-2 


J 


10-^ 


] 




10 


-j 


1 


-j 


10-' 


-j 


10-2 


-j 


10-^ 


■ 




10 




1 


n 


10-' 




10-2 


n 


10"^ 




10 


n 


1 




10-' 




10-2 




10-^ 





neutron PDR 
excitations « 



r 

neutron PDR 
excitations 
V. 



V 












- — T— 




1 


a] 


1 


1 



neutron PDR 
excitations 



1 ' r 

proton PDR 
excitations 



GDR 



GDR 



GDR 



GDR 



Lb 



4 



T ' r 
GDR 



10 



12 



14 



E [MeV] 



16 



132 



Sn 



124 



Sn 



'I' i' I 



122 



Sn 



116 



Sn 



100 



Sn 



18 



] 



20 



FIG. 7: (Color online) QRPA results for the PDR and GDR strength distributions in 

100,116,122,124,132Sn igotopes. 



The good overall agreement between the calculations and the experiment 36(] for the total 
PDR strength and centroid energy in ^^^Sn indicates that in this particular case the main 
PDR properties could already be determined on the level of the one-phonon approximation. 
A similar conclusion was drawn in our previous QPM calculations for several Sn isotopes 
with A=120-130 Gj. An important observation is that the low-energy tail of the GDR can 
give a strong contribution to the dipole strength around the particle threshold. This effect 
appears because the GDR states may overlap with the PDR region and can be fragmented 
due to coupling to multi-phonon states. The effect becomes increasingly important in nuclei 
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where the neutron threshold is higher, hence approaching the GDR region, as in the hghter 
Sn isotopes or ^^^Sn, where the PDR strength is situated very close to the neutron threshold. 
Such a situation demands much larger model spaces. 

VI. OTHER MODEL CALCULATIONS AND EXPERIMENTAL DATA 
A. PDR Models 

Overall, the present Sn results are in, at least qualitative, a gree ment with the theoret- 
ical PDR investigations by Density Functional Theory (DFT) [241, relativistic RPA j4o| . 



relativistic QRPA 







42l |. extended theory of finite Fermi systems J43l| and QRPA-PC 
(Quasiparticle Random Phase Approximation plus Phonon Couplin g') |44l| . They confirm 
the conclusions drawn from our former QPM calculations in ^osp]^ j^j^gj ^j^g N=82 

case, studied recently in Sj. All these different approaches confirm the PDR mode as a 
universal low-energy dipole mode of a character generic for isospin asymmetric nuclei. 



As an example we cite the studies of ref. 



2J] investigating low-lying dipole states in ^^Ca 



in a density functional theory approach. Similar to the much heavier nuclei considered here 
the PDR is predicted to be located in the energy range 5-10 MeV. Also in that nucleus, the 
centroid energy of the PDR strength is found to decrease with the number of the neutrons, 
while the integrated PDR strength (below the neutron particle emission threshold] increases. 
These results agree with the present calculations and our findings in refs. jsl, Q, [s] in the 
Sn isotopes. A common result of all model calculations discussed here is a clear connection 
between the existence of low-energy dipole strength and the presence of isospin assymetry 
or nuclear skin in the investigated nuclei. However, we emphasize, that arguments based on 
the energy alone are likely to be insufficient for a unambiguous identification of the dipole 
states as belonging to the PDR. To our understanding, as an important conclusion from the 
analysis of the transition densities, the PDR strength is attached only to the states located 
below the neutron particle emission threshold. Hence, the centroid of the PDR energy has a 
tendency to be closely connected with the one-neutron separation energy. At higher energies, 
the dipole spectrum merges rapidly into the low energy tail of the GDR and the transitions 
lose their characteristic PDR features. 

A controversial question is the degree of the collectivity of the PDR transitions. This 
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issue has been discussed by several authors [2J, |40|, |4l|, 14^, |45| . In the non-relativistic models 
like ours 3, S [3] and QRPA-PC for example, the PDR is referred to as the excitation of 
two-quasiparticles states. In ^^^Sn the relativistic QRPA 4l| predicts a collective neutron 
state at E* = 8.6 MeV that has been related to the PDR excitation. This state contains 
particle-hole configurations accounting for transitions into continuum states. The collectivity 
of such excitations we found to dependent strongly on the choice of the spin-orbit potential 

Dv shifting the 



46j . otherwise 



affecting the energy gap between the bound hole and unbound particle region 
continuum states. Our standard choice for the spin-orbit potential strength 
describing the spectra reasonably well, disfavors such admixtures. 

Experimental data for low-energy dipole states below the particle emission threshold are 
available for a number of Sn isotopes, ^^^Sn and ^^^Sn jsG] and recent measurements in ^^^Sn 



37|, and for several N=82 isotones Altogether, our calculations describe these data quite 



satisfactory. 



B. Dipole Response in I30,i32gjj 

We pay special attention to the region around ^^^Sn, because of the expected closure of 
the N=82 neutron shell as indicated e.g. by the energy of the first 2+ state. The HFB 
calculations predict a double shell closure for protons and neutrons, respectively. On the 
other hand, the HFB calculations show, that the N=82 neutron shell closure depends to 
some extent on the balance between spin-orbit splitting and the effective pairing strength. 

The LAND-FRS collaboration at GSI has recently measured in a pioneering Coulomb 
dissociation experiment the dipole response above neutron threshold in i30,i32gj^ These 
measurements are providing the first data on the dipole response in the these exotic nuclei. 
However, we have to be aware, that any dipole strength below the particle emission threshold 
- if existing - cannot be accessed by this type of measurement. Besides the GDR a prominent 
feature of the data is a resonance-like structure around E* ~10 MeV exhausting a few percent 

n 

of the EWSR in ^^°'^^^Sn nuclei. In [5| this part of the response function was interpreted 
as a PDR. We compare our calculated integrated dipole photoabsorption cross sections a in 
the Sn isotopes to the LAND-FRS data fl) in Table [ml 

A different conclusion is obtained by analyzing our QRPA wave functions and the dipole 
transition densities. In ^'^^'^'^^Sn, we find dipole excitations, carrying the characteristic fea- 
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TABLE III: Dipole response in ^^'''^^■^Sn. Calculated energies and integrated cross sections (columns 
denoted by QPM) in one-phonon approximation are compared with recent measurements (columns 
denoted by Exp.) [5] of PDR and GDR in Sn isotopes. The calculated integrated PDR and low- 
energy GDR cross sections are denoted by J a^^^ and / c^Jr^ > respectively. The total photoab- 
sorption cross section up to 20 MeV is denoted by J (jGDR 

Nucl. PDR {E)^^^ fa^^^ E^^J f a^^^ E^E^ f a^E^ E^^^^ E^ff^ f a^^^ J a^^^ 
iregTon) [^eV] [mb MeV] [MeV] [mb MeV] [MeV] [mb MeV] [MeV] [MeV] [mb MeV] 
QPM QPM QPM Exp. Exp. QPM QPM Exp. QPM Exp. QPM 
i30Sn 0-7.4 5.8 8.2 10.1(7) 130(55) 8-11 137.3 15.9(5) 16. 1930(300)* 1616 

i32Sn 0-8 7.1 10.4 9.8(7) 75(57) 8-11 97.6 16.1(7) 16.1 1670(420)* 1518 



*Data of ref.^a] integrated up to 20 MeV 



48|. 



tures of PDR transitions, below the neutron particle emission threshold, as indicated in 
the first three columns of Table IIIII and in FigH In the energy domain E'* = 8 — 12 MeV, 
assigned in 5|] as PDR region, we obtain in both nuclei another concentration of El strength 
(see also Fig. [7]). However, because the transition densities show the GDR-type behavior, we 
consider this part of the dipole response as the low-energy tail {LET) of the GDR. Although 
the LET evolves in close relation to the neutron excess, but it does not seem to be related 
to excitations of the neutron skin. 

In fact, there is a simple proportionality between the dipole photoabsorption cross section, 
integrated over an energy interval around E*., and reduced transition strength, j dEa^ ~ 
E*B{E1, E*), up to a numerical factor [34]. Exploiting this relation, we have calculated the 
integrated dipole photoabsorption cross sections in the LET regions of ^^°'^^^Sn. In Tab JIIII 
it is seen that the theoretical results agree rather well with those determined experimentally 
in jsl. Within the experimental error bars, also the full strengths, including excitations up 
to 20 MeV, are reasonably well described. 

For the purpose of a realistic description of the measured spectra, we have applied a 
slightly different numerical method by solving the QRPA Dyson equation similar to the 



approach used in 



471 ] allowing to take into account explicitly the continuum decay width T^ 



of the states above particle threshold, ranging from a few keV up to about hundred keV. 
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FIG. 9: Electromagnetic QRPA dipole response function and LAND-FRS data [5(1 for ^'^^Sn. The 
QRPA results (left), where obtained by solving the Dyson equation and include the decay widths 
from particle emission. In the right panel the theoretical response function has been folded with 



the experimental acceptance filter 



48( 1 (dashed line) and is compared to the data (symbols). 



Still, a comparison to the LAND-FRS 



the experimental acceptance filters 



spectra is only possible after folding the theory with 



481 ]. The results of such calculations is shown in Fig. [9] 



and FiglH] with a quite remarkable agreement to the data. 



VII. SUMMARY AND CONCLUSIONS 

Low-energy dipole excitations in ^^^^^^"^871 isotopes were studied by a theoretical approach 
based on HFB and QPM theory. From our calculations in ^^^Sn ^^^Sn we obtained low- 
energy dipole strength in the energy region below 8 MeV, close but below the neutron 
emission threshold. These states are of a special character. Their structure is dominated 
by neutron components and their transition strength is directly related to the presence of a 
neutron skin. Their generic character is further confirmed by the shape and structure of the 
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FIG. 10: Electromagnetic QRPA dipole response function and LAND-FRS data [5] for ^^^Sn. As 
in Fig. [9] the QRPA results (left) include the decay widths from particle emission. In the right 
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panel the theoretical response function has been folded with the experimental acceptance filter 
(dashed line) and is compared to the data (symbols). 

related transition densities, showing that these pygmy dipole resonances (PDR) are clearly 
distinguishable from the conventional GDR mode. Our calculations show a rather abrupt 
transition from PDR- to the GDR-type excitations, typically occurring at energies slightly 
above the particle threshold. An important finding in our calculations is that an accurate 
description of the PDR part of the dipole spectrum requires a single particle spectrum 
corresponding to a total effective mass ^ = 1. From that observation we conclude that non- 
localities and dynamical effects form core polarizations are important for a proper description 
of the PDR spectrum. Pure HF and HFB models, whether non-relativistic or relativistic, 
typically use effective masses considerably less then unity. Hence, such approaches might 
miss important effects. 

In the most proton-rich exotic nuclei ^00-104 g^^ ^j^g lowest dipole states are almost pure 
proton excitations. They are related to oscillations of weakly bound protons, indicating a 
proton PDR. The interesting point is, that we found these states in heavy nuclei with N 
slightly larger, or equal to Z. We suggest, that the effect is due to Coulomb repulsion, that 
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pushes out weakly bound protons orbitals to the nuclear surface. Coulomb effects also induce 
a considerable amount of isospin breaking at the level of single-particle wave functions. 

Since similar observation have been made in the nearby N = 82 isotonic nuclei, we may 
conclude, that the features discussed here are indicating a new universal mode of excitation. 
It is worthwhile to extent the investigations also in other mass regions. Promising candidates 
are the Ni and Ca isotopes, but also the light mass region, where a mixing between halo and 
skin degrees of freedom can be expected, which may lead to still other modes of excitations. 
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APPENDIX A: THE REARRANGEMENT POTENTIALS 

Once the proton and neutron self-energies Sp ,j(p), respectively, are known the rearrange- 
ments parts arc determined and properly subtracted by exploiting relations found in infinite 
nuclear matter. In symmetric nuclear matter with pp — p„ we find for the isoscalar self- 
energy Eo = (S„ -I- Sp)/2 the relation 

So(p) = \§pPMPl (Al) 

which we integrate to give 

C/o(p) = - rdp'Eoip') , (A2) 
pJo 

providing us with Uo{p) — {Un{p) + Up{p))/2. In pure neutron matter we have p — pz — Pn 
and 

Un{p) = - r dp'T.n{p') . (A3) 



PJo 

This allows to determine 



Up{p) = 2U^{p)-Un{p) . (A4) 
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For a finite nucleus the densities are given parametrically as functions of the radius r. Hence, 
we can replace the integrations over density by radial integrals 



where p{r) is the density calculated self-consistently according to eq. 13.31 with wave func- 
tions from the effective potential ^^(r). Obviously, the above equation is applicable to any 
potential given as a function of the radius. Hence, by means of these defolding relations we 
are able to calculate B{A) for arbitrary phenomenological single-particle potentials, which 
otherwise we could not. 
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